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We demonstrate the remnant presence of initial correlations in the steady-state electrical current 
flowing between low-dimensional interacting leads. The leads are described as Luttinger liquids 
and electrons can tunnel via a quantum point-contact. We derive an analytic result for the time- 
dependent current and show that ground-state correlations have a large impact on the relaxation 
and long-time behavior. In particular, the I-V characteristic cannot be reproduced by quenching 
the interaction in time. We further present a universal formula of the steady-state current js for an 
arbitrary sequence of interaction quenches. It is established that js is history dependent provided 
that the switching process is non-smooth. 
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Non-equilibrium phenomena in open nanoscale sys- 
tems offer a formidable challenge to modern science [1]. 
Controlling the electron dynamics of a molecular de- 
vice is the ultimate goal of nanoelectronics and quantum 
computation[2]; its microscopic description a problem at 
the forefront of statistical quantum physics[3]. Resorting 
to approximate methods is inevitable to progress. 

Standard many-body techniques consider an initial 
state with no interaction and with no contact between 
the system and the baths (leads from now on) , and then 
switch them on in time [4-6]. In fact, it is plausible to 
believe that starting from the true interacting and con- 
tacted state the long-time results would not change. To 
what extent, however, such belief is actually the truth? 
This question is of both practical and fundamental inter- 
est. It has been shown by us [7] and others [8] that for 
non-interacting electrons the initial contact plays no role 
at the steady-state [9] (theorem of equivalence) . Allow- 
ing for interactions in the system only (non-interacting 
leads) Myohanen et al. found that steady-state quanti- 
ties are not sensitive to initial correlations either [10]. It is 
the purpose of this Letter to show that interacting leads 
change dramatically the picture: the switching process 
can indeed have a large impact on the relaxation and the 
steady-state behavior. 

We consider two one-dimensional interacting leads de- 
scribed as Luttinger Liquids (LL)[11], see Fig. 1. It is 
known that a LL does not relax to the ground state after a 
sudden quench of the interaction [12-14] (thermalization 
breakdown) . The implications of such important result in 
the context of time-dependent (TD) transport are totally 
unknown and will be here explored for the first time. We 
compare the dynamics of initially (a) contacted (t]t = 1) 
versus uncontacted (r/r = 0) and (b) interacting (rji = 1) 
versus non-interacting (r/i = 0) LL when driven out of 
equilibrium by an external bias. Our main findings are 



FIG. 1. Sketch of the device. Two interacting leads hosting 
L and R movers are connected at x = via a weak link. A 
bias voltage Vl — Vr can be applied between the leads. 



that in (a) the system relaxes towards the same steady- 
state although with a different power law decay. In (b) 
the sudden quench of the interaction when rjj = al- 
ters the steady- current js as well. This remains true for 
an arbitrary sequence of interaction quenches. We are 
able to write js as an explicit functional of the switching 
process and to establish that js is history dependent for 
non-smooth switchings. 

The equilibrium Hamiltonian for the system of Fig. 1 
reads 



Ho = H R + H L + rjiHi + ri T H T . 



(1) 



The one-body part of the left (L) and right (R) leads 
is H R/L = ^fivpj^dx tp] i/L (x)d x if! R/L (x), where the 
fermion field ipR./L describes right/left moving electrons 
in lead R/L with Fermi velocity vf (chiral leads). We 
take a density-density interaction of the form Hj = 
U-oo dx l 2 d2Pn(x)pL{x) + 9i{p R {.x) + pl(x))}, where 

p R ji 1 = : WniL^R/h '■ is (in standard notation) the 
fermionic density operator relative to the Fermi sea, and 
52/4 are the forward scattering couplings, corresponding 
to inter/intra lead interactions respectively. The two chi- 
ral liquids are linked at x = via the tunneling term 
Ht = A-0}j(O)'(/'l(O)-|-H.c., which docs not commute with 



the total number of electrons Nr/l of each lead. 

If a bias V = Vj, — Vr is applied at, say, time £ = 0, a 
finite current j(t) starts flowing across the link. The cur- 
rent operator (in atomic units) J = dN^/dt = —dN^/dt 
reads J = iAV'|j(0)V , l(0) +H.c. At zero temperature the 
current j(t) is the TD average of J over the ground state 
|*o) of H , i.e., 



j(t) = (<H \J Hl (t)\^ ), 



(2) 



where J^ (£) is the J operator in the Heisenberg repre- 
sentation with respect to the interacting, contacted and 
biased Hamiltonian Hi = Hl + Hr + Hj + Ht + Hy, 
Hy = VrNr + VlNl. Note that the factors t]i,t]t refer 
to times t < and different values r\i,i]T = 0, 1 yields 
different Ho and hence different initial states |\&o)- At 
positive times the Hamiltonian is the same in all cases. 

The exact non-interacting solution. We start our anal- 
ysis by calculating j(t) when tjt — (initially uncon- 
tacted) and gi — 54 = (always non-interacting). 
In terms of the Fourier transform tpk R/L of the orig- 
inal fcrmion fields, the current operator reads J = 
(iX/a) X^fcfe' Yk ipk' + H.c., with a the usual short- 
distance cutoff. Its expectation value is then j(t) = 
AIm E Q /v r p Q W/p [ r p Q W]* where the sum runs 



over a — R,L, /, 



R/L 



f(±VFp) is the Fermi function of 



lead R/L and Tf(t) = -to/ f (^ol^e - "* 1 *^!*!)) 
is the sum of the probability amplitudes (retarded 
Green's functions) for the transition pp — > k a . From the 
Dyson equation it is straightforward to find T" a (t) = 
_ ie i(av F p+v a )t J (i + c 2) and r« a (t) = -icT« a (t), with 
c = \/(2vf), and hence 



m 



2c 2 



7r(l 



„2\2 



V. 



(3) 



The current is discontinuous in time; the steady-state 
value is reached instantaneously. This is due to the un- 
bound (relativistic) energy spectrum [5] and the lack of 
interactions, as discussed in detail in Rcf. 15. As we 
shall sec, when Hi ^ the transient regime is more com- 
plex. 

Current to lowest order in A. Thc problem does not 
have an exact solution when both Hi and Ht are present. 
Below, we calculate j(t) to lowest order in A. In general, 
perturbative treatments in the tunneling amplitude are a 
delicate issue[16]. In our case j(t) has a Taylor expansion 
with convergence radius A < 2vp for Hi = 0, see Eq. 
(3). We, therefore, expect a finite convergence radius at 
least for small interaction strengths. Let the unperturbed 
Hamiltonian be Ho = Hr + Hl + r/iHi in equilibrium 
(t < 0) and Hi = Hr + Hl + Hi + Hy at positive times. 
At zero temperature and to lowest order in A 



j'(*)=*<*o|/ ds H TAi {s),J Al (t) -vr dr 

Jo L J Jo 



with l^o) the ground state of Ho- The first term in the 
r.h.s. is the standard Kubo formula. Such term alone 
describes the transient response when the contacts are 
switched on at t = (t]t = 0). If, however, the equi- 
librium system is already contacted (r/T — 1) we must 
account for a correction; this is the physical content of 
the second term [17]. At any finite time initial correla- 
tion effects are visible in both terms due to the ground 
state dependence on r/i. When t — > 00 only thc Kubo 
term survives, which yields the steady-current js- The 
dependence of js on the ground state (rji = 0, 1) will be 
addressed below. 

The averages in Eq. (4) can be explicitly calculated 
by resorting to the bosonization method[ll]. We intro- 
duce the scalar fields (ft and 9 from pn{x) + pl{x) = 
-^d x (ft{x) and ift R/L (x) = ^^e'v^wT'M], with 
k r/l the anticommuting Klein factors. The scalar 
fields obey the commutation relation [(ft(x),8(x')] — 
zsgn(a; — x')/2. In terms of (ft and 9 the Hamilto- 
nian H = Hr + Hl + Hi is a simple quadratic form 
H =JJ^dx[K-Hd x( ft(x)) 2 + Kd x 6(x)) 2 }, with v = 

y // (2irvF + gi) 2 — g'i)/^ 71 the renormalizcd velocity and 
K = yj (2-kvf + <?4 — g2)/(2irvF + .94 + 32) a parameter 
which measures the interaction strength. Note that 
< K < 1 for repulsive interactions; K = 1 corresponds 
to the noninteracting case while small values of K indi- 
cate a strongly correlated regime. 

By employing the gauge transformation[18] ipL,R ~^ 
tpL,Re tVL/Rt the problem of evaluating Eq. (4) is re- 
duced to the calculation of different bosonic vacuum 
averages [11]. After some tedious algebra one finds 



j(t)=^Rc[r lT A vl (t) + B VI (t)} 



(5) 



where 



A (t) = sm(Vt) / dr 7 2 (t + It) , 
Jo 

B (t) =i ds sm[V(s - t)] j 2K (s - t) 



x | 7 (s-£)| (1 - K)2 
for r/i = and 



7 2 ( S + t) 



7(2t)7(2a) 



l-K^ 



(6) 



Ai{t) =sin(V£) / drY (t + ir) 
Jo 

Bi(t) = i I ds sm[V(s - t)} j 2K (s - t) 
Jo 



(7) 



for rji = 1, and where 7(2) = a/(a — ivz) and £ = 
X 2 /{-Ka) 2 . In all cases (rn,r]T = 0,1) j(t) is an odd 
function of V, as it should be. We also notice that for 
noninteracting systems (K = 1) we recover the expected 
result Ai = Ao and Bi —Bo- In this case the function 
^RefB^o] coincides with the current in Eq. (3) to lowest 
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FIG. 2. Transient currents Jti (t) (solid) and JTo(t) (dashed) 
for V = 10" 2 , K = 0.88 (upper panel), and K = 0.75 (lower 
panel). In the long-time limit they reach the same steady- 
state value. Current is in units of £a/v, V is in units of v/a 
and t is in units of W^a/v. 



FIG. 3. Transient currents jn(t) (solid) and jio(t) (dashed) 
for V = 10" 2 , K = 0.75 (upper panel) and K = 0.88 (lower 
panel). In the long-time limit they reach different steady- 
states. Same units as in Fig. (2). 



order in A. We can now provide a quantitative analy- 
sis of the TD current response for different preparative 
configurations. 

Contacted versus uncontacted ground state. We con- 
sider an initially contacted (j]t = 1) and uncontacted 
(tjt = 0) correlated ground state (rjj = 1) and compare 
the corresponding TD currents jti = £Re[Ai + B{\ and 
jto = £Re[.Bi]. The current JTo(t) has been recently 
computed in Ref. [19]. In the long time limit it returns 
the well known steady-state result 



j s (P)=sm(nK) K ([3)sgii(V)\Vf 



(8) 



with k((3) = -^(a/w)' 3 + 1 r(-/3) sin(/37r/2) and the ex- 
ponent /? = 2K — 1, obtained long ago by Kane and 
Fisher [20]. Since A\(t — > oo) = 0, jti approaches the 
same steady state. Note that the small bias limit is ill- 
defined for K < 1/2 due to the break down of the per- 
turbative expansion in powers of A[18, 21]. Even though 
JTo{t — > oo) = JTi(t — > oo) the relaxation is different 
in the two cases, see Fig. 2. The function JTo(t) ap- 
proaches the asymptotic limit with transient oscillations 
of frequency V and damping envelope proportional to 
t~ 2K [19]. The more physical current jti, instead, de- 
cays much slower. The integral in A\ (t) can be calculated 
analytically and yields 



JTi(t)-JTo(t) =£a 



2 ^-sin(Vi) cos[(2.K~ — 1) arctan(v£/a)] 



2v(2K-l)(a 2 +v 2 t 2 ) K - 1 / 2 



(9) 

which for long times decays as t 1 2K . (Equation (9) pro- 
vides an independent, TD evidence that the perturba- 
tive treatment breaks down for K < 1/2.) Thus, an ini- 
tially contacted state changes the power-law decay from 
~ t~ 2K to the slower ~ t 1 ^ 2K . The amplitude of the 
transient oscillations is also significantly different, due to 
the factor (2K- 1)- 1 in Eq. (9). For K = 0.75, j T1 oscil- 
lates with an amplitude about 10 times larger than that 



of Jto, see Fig. 2. The magnification of the oscillations 
was unexpected since for jti we only switch a bias while 
f° r Jto a l so the contacts. This effect is not an artifact 
of the perturbative treatment: to support the validity of 
our results we checked that for t]i = r\T = 1 and zero 
bias the density matrix p(t) = (^owrh (*)V'i,Hi(i)|^ r o) 
does not evolve in time to first order in X (this is obvi- 
ous for the exact density matrix). The constant value 
p(i) = p(0) is the result of a subtle cancellation of TD 
functions similar to Ai(t) and B\{t). 

Correlated versus uncorrelated ground state. Next we 
consider the effects of correlations in the ground state. 
We take tjt — 1 and compare the TD currents jn and 
jio resulting from Eq. (5) when r\i = 1 and rjj = respec- 
tively. Note that jn = jti (already calculated above). 
The current jio = £Re[^4o + -Bo] is the response to a sud- 
den bias switching and interaction quench; at t > the 
electrons start tunneling from L to R and at the same 
time forming interacting quasiparticles. The interaction 
quench has a dramatic impact on the transport prop- 
erties, both in the transient and steady-state regimes. 
From Fig. 3 we clearly see that the relaxation behavior 
is different. The damping envelope of jioit) is propor- 
tional to t~ K _1 as opposed to t 1 ~ 2K of jn(t)- Notice 
that the exponent ~K 2 — 1 < for all K (first-order 
perturbation theory in A is meaningfull for all K). 

In the long-time limit we find the intriguing result that 
jio(t — > oo) is exactly given by Eq. (8) with exponent 
(3 = K 2 , thus suggesting that the structure of the for- 
mula (8) is universal. Below we will prove that this is 
indeed the case and that /3 is an elegant functional of the 
switching process. For now, we observe that ground state 
correlations are not reproducible by quenching the inter- 
action. The system remembers them forever and steady- 
state quantities are inevitably affected. This behavior 
is reminiscent of the thermalization breakdown enlight- 
ened by Cazalilla[12] and others[13, 14]. Here, however, 
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FIG. 4. Transient currents with r/i = 1, r\T = (solid), 
and with ni = tjt ~ for the quench 1 — ^ if (dashed) and 
the quench sequence 1 — > LhS. — >■ if (dotted-dashed). Here 
if = 0.75 and 1/ = 10 -2 and the second quench occurs at 
£i = 1; same units as in Fig. (2). 



we are neither in equilibrium nor close to it (the bias 
is treated to all orders). The non-equilibrium exponents 
/3 = 2K — 1 and (3 = A 2 refer to current-carrying states 
as obtained from the full TD Schrodingcr equation with 
different initial states. 

History dependence. We now address the question 
whether or not the physical steady-current js(2K — 1) of 
Eq. (8) is reproducible by more sophisticated switching 
processes of the interaction like, e.g., an adiabatic switch- 
ing. Preliminary insight can be gained by calculating j(t) 
for a double quench: we first quench an interaction with 
K\ = (I + K)/2, let the system evolve, and then change 
suddenly K\ — > Ki = K. The current is calculated along 
the same line of reasoning of Eq. (4) , although the formu- 
las become considerably more cumbersome. In Fig. 4 we 
compare the TD currents for initially uncontacted leads 
resulting from an interaction A (solid), a single quench 
1 — > K (dashed), and the aformentioned double quench 
(dotted-dashed). We clearly see that in the latter case 
the steady-current is larger than js(K 2 ) (single-quench) 
and gets closer to js(2A — 1). Strikingly, the double- 
quench steady-current is again given by js(P) of Eq. (8) 
with p = |(1 + Af)(l + (f^) 2 ) - 1. This value of (3 
depends only on the A-sequence and is independent of 
the quenching times. We have been able to extend the 
above solution to systems initially interacting with K$ 
and then subject to an arbitrary sequence of quenches 
Kq — >• K\ — > . . . — > Km = A. We found the remarkable 
result that the formula (8) is universal, with the sequence 
dependent j3 given by 



JV— l 



Wn] ^n 



1 



A„ +] 
~K^ 



- 1. 



(10) 



This formula yields the correct values of f3 for the sin- 
gle and double quench discussed above. Note that for 
a sequence of increasing interactions K n +i < K n it 
holds j3 > 2K — 1 with the equality valid only for 
K = K 1 = ... = K N = K. 

We now show that the special value /3 = 2K — 1 is also 
reproducible by an arbitrary (not necessarily adiabatic) 



continuous (N — ¥ oo) sequential quenching. In this limit 
the variable x n = n/N becomes a continuous variable 
and we can think of the K n as the values taken by a dif- 
ferentiable function K(x) in x — x n , with K(0) = K and 
K(l) = K . Then, the quantity j3 becomes a functional 
of K (x) that we now work out explicitly. Approximating 
K(x n+1 ) = K(x n + jj) ~ K(x n ) + jjK'(x n ) and taking 
the logarithm of Eq. (10) we can write 



log 



( WW 



V 2#(0) 



log 1 



N-l 

= lim y 



1 K'{x n ) 

N K(x n ) 
K{1) 



/0 K(x) °K(oy 

from which it follows the history independent result 
P[K(x)] = 2 A - 1. 



(11) 



(12) 



The above result can easily be generalized to discontin- 
uous switching functions K(x) for which, instead, the 
exponent j3 is history dependent. 

Conclusions. In conclusion we studied the role of dif- 
ferent preparative configurations in TD quantum trans- 
port between LLs. By using bosonization methods we 
showed that a sudden switching of the contacts does 
not change the steady-state but alters significantly the 
transient behavior, changing the damping envelope from 
~ t~ 2K to ~ t 1 ~ 2K and magnifying the amplitude of the 
oscillations. The effects of a sudden interaction quench 
is even more striking. Besides a different power law de- 



cay (~ t 



1-2K 



versus 



t 



-K z 



" x damping envelope) the 



steady-current is also different; the TV characteristic 
js oc y 3 changes from /3 = 2 A — 1 to j3 — A 2 . More 
generally we proved that for a sequence of interaction 
quenches the steady-current is a universal function of the 
exponent /? which, in turn, is a functional of the switch- 
ing process. It is only for smooth switchings that (3 is 
history independent and equals the value 2A — 1 of the 
initially interacting LL. The explicit j3 functional derived 
in this Letter establishes the existence of intriguing mem- 
ory effects that point to a complex entanglement between 
equilibrium and non-equilibrium correlations in strongly 
confined systems. 
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